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ABSTRACT 

We have used the Xco conversion factor, Local Thermal Equilibrium and Large Velocity 
Gradient approximation to parametrize the cold and warm phase of the interstellar medium 
from five different low transitions of the CO molecule in the central 21" (kpc) region of 
NGC 1097. We have applied a one-component model and derived a typical kinetic tempera- 
ture of about 33 K, a molecular Hydrogen density of 4.9 x 10 3 M Q pc~ 3 and a CO column 
density of 1.2 x 10~ 2 M pc~ 2 . A two-component model results in 85% cold-to-total gas 
fraction in the presence of a 90 K warm counterpart. Furthermore, we "resolve" the spatially 
unresolved single dish observations by selecting velocity channels that in an interferometric 
velocity map correspond to specific regions. We have selected five such regions and found 
that the physical properties in these regions are comparable to those derived from the full 
line profile. This implies that the central kpc of NGC 1097 is rather homogeneous in nature, 
and, although the regions are not uniquely located within the ring, the star formation along 
the ring is homogeneously distributed (in agreement with recent Herschel observations). We 
have further revised the mass inflow rate onto the Supermassive Black Hole in this prototype 
LINER/Syl galaxy and found that, accounting for the total interstellar medium and applying 
a careful contribution of the disc thickness and corresponding stability criterion, increases 
the previous estimations by a factor 10. Finally we have calculated the Xco conversion fac- 
tor for the centre of NGC 1097 using an independent estimation of the surface density to the 
CO emission, and obtained X C o = (2.8 + 0.5) x 10 2Q cm -2 (Kkms -1 ) -1 at radius 10.5" and 
^co - (5-0 ± 0.5) x 10 20 crrT 2 (Kkms -1 ) -1 at radius 7.5". With the approach and analysis de- 
scribed in this paper we have demonstrated that important physical properties can be derived 
to a resolution beyond the single dish resolution element, however, caution is necessary for 
interpreting the results. 

Key words: Galaxies: individual: NGC 1097, Galaxies: interstellar medium, Galaxies: fun- 
damental parameters 



1 INTRODUCTION 

The connection between star formation in centres of galaxies and 
the onset of nuclear activity involve processes still not completely 
understood. In many cases, the majority of star formation in galax- 
ies is located in circumnuclear rings, and the most extended mech- 
anism for fuelling them and generate densities high enough to form 
stars, is believed to be non-axisymmetric perturbations, such as bars 
or tidal interactions. In the presence of a bar, over-densities of gas 
may occur in resonances, like Inner Lindblad Resonance (ILR) for 
a inner ring (e.g. |Combes & Gerin|1985||Athanassoula|1992||Wada| 
\l99ty . Once the gas is accumulated in the ring, there are two main 
hypothetical scenarios for the formation of the stars in the ring: (1) 
star formation produced by gravitational instabilities in different gi- 
ant clouds placed randomly along the ring (Elmegreen 1994), and 
(2) quasi-instantaneous star formation in two regions with over den- 



sities, usually the connection points of the dust lanes of the bar with 
the ring, followed by migration of the star forming knots along the 
rotation velocity of the disc, producing an age gradient along the 
ring ( |Boker et al.|2008) . 

The study of star formation locations is typically done observ- 
ing Ha or radio continuum emission. In addition, it is possible to 
use direct measurements of gas column density (E gas = Shi + £h 2 ), 
to which the star formation rate per unit area X S fr is empirically 
related by the Schmidt law ( |Schmidt| 1 959| |Kennicutt| 1 989) . One 
aim of this paper is to measure the molecular gas column density in 
order to characterize the star formation along the 1 kpc ring of the 
nearby prototype LINER/Syl barred galaxy NGC 1097, estimated 
through CO observations. 

The fuelling of nuclear and/or circumnuclear activity in barred 
galaxies is driven by the bar itself due to loss of angular momentum 
along the dust lanes (e.g. [Athanassoula 1992 ), however, the fate 
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of the inflow mass after passing the ILR is still unsolved. Several 
mechanisms have been inferred to explain the fate of the inflowing 
gas. One proposed mechanism is the coupling of nested bars inside 
the ILR of primary bars, transporting gas and stars towards the cen- 
tre (e.g. |Shlosman et al.|1989]|Englmaier & Shlosman|2004) . This 
mechanism is not completely established since there are other stud- 
ies which show that nested bars are not allow to carry material to 
the inner few parsec of galaxies (Maciejewski et al. 2002). Another 
proposed mechanism is inflow along nuclear spiral arms, which are 
connected with the dust lanes of the primary bar (e.g. Englmaier 



|& Shlosman1[2000l |Maciejewski et al.||2002| |Maciejewski||2004 1 
These arms have a very low density contrast compared with the am- 
bient inter-arm regions, only about the 10% (Englmaier & Shlos- 
|man|20 00l, and present a wide variety of observed morphologies. 
Morphologically, nuclear spirals are divided into two main cate- 
gories, chaotic spirals, and symmetric or grand design spirals, for 
which different formation mechanisms have been proposed (e.g., 



Elmegreen et al. 1998, Englmaier & Shlosman 2000, Maciejewski 



et al.|2 002). In the presence of a central super massive black hole 
(SMBH), spiral shocks can extend all the way to its vicinity and 
generate gas inflow compatible to the accretion rates observed in 
local active galactic nuclei, (e.g. morphologically by |Martini et al.| 
[2003] and kinematically by |Fathi et al.|2006fr . 

NGC 1097 is an excellent candidate to study the inflow pro- 
cesses driven by the bar, however, the lack of direct measurements 
of the physical conditions of the interstellar gas in the circum- 



nuclear regions has complicated matters. Optical Ivan de Ven & 



Fathi|2010| l; Near infra-red ( Davies et al. 2009 and Kotilain en et al 



2000), radio continuum (Beck et al. 2005 ), CO transitions with ra- 



dio interferometers and single dish (Kohno et al. 2003, Petitpas 



& Wilson 2003 Hsieh et al. 2008) and far infra-red observations 
(Beirao et al. 2010 Sandstrom et al. 2010) have been performed 
for its 1 kpc starburst ring. Moreover the still unclear processes in- 
volved in the star formation in the ring, the estimation of the mass 
accretion rate toward its (1.2 ± 0.2) x 10 s M Q SMBH (jLewis & Er~- 



acleous 2006 ) is to date inconclusive. As a second goal, we shed 
light on the total gas mass in the region inside the ring, and ap- 
ply several methods to better understand the physical state of the 
cold+warm (< 150 K) gas in the circumnuclear kpc of NGC 1097. 

We use archival data combined with our own single dish ob- 
servations of the molecular gas in the circumnuclear regions using 
the Atacama Pathfinder Experiment (APEX) telescope. The data 
are used with two different line ratio modelling methods to obtain 
a new estimation of the physical properties of the gas and the ac- 
cretion mass rate towards the SMBH to compare with the previous 
estimations. Moreover, we derive the molecular gas physical con- 
ditions along regions with comparable line of sight velocities, and 
in this way override the limited spatial resolution of the single dish 
observations. In section [2] we present the observations and data 
reduction. In section [3] we explain all the analysis done for esti- 
mating the line ratios, as for study the physical properties using the 
X C o conversion factor, a local thermal equilibrium and a large ve- 
locity gradient approximation for the inner kpc of NGC 1097. In 
section]?] we explain the analysis made within the ring and discuss 
the result obtained there, in section [5] we discuss the mass inflow 
rate towards the inner parsecs, and in section[6]we present a study 
of the real Xco conversion factor value for this galaxy. Finally, we 
conclude in section|7] 



2 OBSERVATIONS AND DATA REDUCTION 

On the 27th of June of 2008, we observed 12 CO (3-2) and 
B CO (3-2) in one pointing at the centre of NGC 1097 with the 
Atacama Pathfinder Experiment (APEX) on Llano de Chajnantor, 
Chile, using the APEX-2 receiver. The Half-Power Beam Width 
(HPBW) is 18" at 345.8 GHz and 18.8" at 330 GHz. In addition, 
we mapped 12 CO (3-2) along the primary bar of NGC 1097 on the 
7th of July of 2009, using the same set up. We mapped the bar with 
an array of 9x3 points. We scaled the antenna temperature adopting 
the value for rj MB presented in the APEX web page of 0.73. Dur- 
ing the reduction of the spectra, we subtracted a first order baseline 
from the profiles. Lines are shown in Fig.^ 

Our data have been combined with published single dish and 
interferometric data from the centre of this galaxy: 12 CO (1-0) 
from |Kohno et al.| < [2003"l > with a HPBW of 15" observed with the 
Nobeyama Radio Observatory, 12 CO (2-1) and 13 CO (2-1) from 
|Petitpas & Wilson| (2003) with a HPBW of 21" observed with 
James Clerk Maxwell Telescope, and l2 CO (1-0) observed with 
Nobeyama Millimetre Array ( |Kohno et al.||2003) . All the auxil- 
iary data were convolved to the common beam size of 21", and 
we discuss 12 CO (1-0) case separately in the following subsection. 
For 13 CO (3-2) transition we used to estimate the ratio between 
12 CO (3-2) and 13 CO (3-2), both with a beam size of 18", we as- 
sume that the ratio changes only marginally between 18" and 21" 
beam sizes. We further observed the 13 CO (3-2) in additional nine 
positions around the central 18.8", with the APEX telescope, 8th of 
October of 2010 and confirm the validity of the bottom right panel 
of Fig.[T]as representative of a 21" observations. 

2.1 12 CO (1-0) 

We combined the 12 CO (1-0) single dish line with the 12 CO (1-0) 
interferometric data from |Kohno e t al. ( 2003 1 to estimate the miss- 
ing flux. The interferometric data were thus convolved to the res- 
olution of the single dish observations and a spectral line from the 
central beam of the galaxy was extracted, resulting in a missing flux 
of 7.7%. Given the relatively marginal missing flux (< 10%), emis- 
sion lines corresponding to 21" beam could be obtained simply by 
co-addition of the interferometric data. 

Furthermore, we improve our estimation of the line using the 
average of two different methods to recover the 21" line. Method 1: 
adding 7.7% to the co-added interferometric data, and Method 2: 
by calculating the missing flux for every channel of the line with 
15" beam, and adding the corresponding missing flux (per channel) 
to the corresponding channel of the generated 21" line. Following 
these tests, for our analysis, we use the average of these two cases, 
for which a flux error can be obtained as the mean of the standard 
deviations of every channel. 



3 THE UNRESOLVED CENTRAL KPC 
3.1 Line-ratio estimation 

All lines were binned to the velocity sampling of 26 km s -1 to 
match the largest sampling available to us. We neglected channels 
with a signal-to-noise (5 IN) smaller than 3, and calculated the ratio 
of the different lines by averaging the ratios between temperatures 
of every channel. The corresponding errors were then the standard 
deviation of this average. We also estimated the ratios using the 
velocity-integrated intensity of every line, and found that they agree 
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Figure 1. The different observations with corresponding transition indicated in each panel at 21" resolution (black curves) as described in section [2] The 
original 15" 12 CO (1-0), 18" 12 CO (3-2) and mapped I3 CO (3-2) lines are over-plotted with red dotted curves. 



Table 1. Spectral line parameters and specifications. 



Line 


Beam size (HPBW) 


Intensity / T mb dV 

Kkms- 1 


Reference 


I2 CO(1-0) 


21/15 


89 ± 9/ 96 ± 2.8 a 


|Kohnq et al.||2003) 


12 CO (2-1) 


21 


110 ±7 


Petitpas & Wilson|{2003J(T2T± 1 K km s^ 1 )" 


12 CO (3-2) 


21/18 


91 ±6/102 ±4 


APEX data b 


13 CO (2-1) 


21 


14 ±2 


Petitpas & Wilson||2003| (14.6 ± 0.6 K km s~') a 


13 CO (3-2) 


18/21 


7 + 3/6 + 3 


APEX data 



a Previous estimation of the intensities by corresponding author. 
b |Petitpas & Wilson|j2003| derived 121.7 ± 0.5 K km s _1 for 21" beam. 



very well with the previous measurement. The results are shown in 
the tablesQ]and[2] 

The errors which contribute to the line measurements are: (i) 
the rms, (ii) the error due to bad pointing, (Hi) the missing flux 
of 12 CO (1-0) and (iv) the uncertainty of the main beam efficiency 
77 MS and that of the Main Beam size 6 M b- Our data reduction and re- 
sampling method provides a good estimation of the rms and miss- 
ing flux errors, and the errors due to beam size and efficiency, to- 
gether, are estimated to be of the order of 10% (at most 10% for 
the main beam efficiency and 0.05% for the main beam size). The 
error due to bad pointing was calculated using the interferometric 
data for 12 CO (1-0) and the error pointing of each telescope (2" 
rms for APEX and 3" rms for the data from Petitp as~& Wilson] 
(2003 1). We studied the difference of the lines produced by varying 
the centre within the error in the pointing, and find that this effect 
corresponds to additional 11% for 2" and 15% for 3" rms. 

3.2 Local Thermal Equilibrium and X C o Conversion Factor 
Analysis 

With the estimated intensities and ratios of the lines, we used two 
different methods to estimate the molecular gas physical properties 
in the inner kpc (21") of NGC 1097, (i) the X C o conversion factor 
and (ii) the Local Thermal Equilibrium (LTE) approximation. Al- 
though, none of these methods are fully realistic for our unresolved 
data, they are instructive to compare with the estimations using a 
Large Velocity Gradient (LVG) approximation. 

Using the Xco conversion factor and the line I2 CO (1-0) we 



estimated the H 2 column density of the inner 21". We used a value 
for Z co of 3 x 10 20 cm ^Kkms 1 ) 1 ( |Scoville et all|l987| and 
|Solomon et al.|1987| also used by |Hsieh et al.|20 08 1 for an illus- 
trative estimation of the molecular gas density. With our estima- 
tion of the integrated intensity of the line, we obtained a column 
density of (2.7 ± 0.3) x 10 22 ctrT 2 . This value is of the same or- 
der as the estimation made by Hsieh et al. 2008 at the radius of 
the starburst ring in this galaxy. However, using the conversion 
factor obtained from outside the Galactic plane Xco = 1.8 x 10 20 
cirr 2 (Kkm s -1 ) -1 {Dame et al.|2001 1 implies a column density of 
(1.6 ± 0.2) x 10 22 cm 2 , and X co = 0.4 x 10 20 cirT 2 (Kkm s" 1 )-' 
I Martin et al. 2010| > gives a column density of (3.6±0.4)x 10 21 cirr 2 . 

We also analysed the emission assuming LTE, making use of 
a "Population Diagram" (Goldsmith & Langer 1999), shown in 
Figure [2] The LTE approximation gives a linear relation between 
ln(Nj/gj) and Ej/k B for a cloud, where Nj is the column density 
of the molecules with the energy state j (Ej) and g,- the statistical 
weight for Ej. The slope of this linear relation, is inversely pro- 
portional to the kinetic temperature characteristic of the cloud. We 
estimated a typical kinetic temperature for 12 CO and 13 CO of about 
10 K. In the case of an optically thick medium, this approxima- 
tion only yields a lower limit for the column density. We derive a 
value that is indeed too low, which is due to the fact that in our LTE 
approximation we do not account for optical thickness. 
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Table 2. Observed line ratios with average errors. 



'■'CO(l-O) 
I2 CO(2-l) 



12 CO(3-2) 
12 CO(2-l) 



12 CO(2-l 



1: CO(3-2) 
I3 CO(3-2) 



Average temperature ratios 
Velocity-integrated line ratios 



0.77 ±0.15 
0.8 ± 0.4 



0.8 ± 0.2 
0.8 ± 0.3 



7.9 



: 1.6 
:4 



18 ± 10 
14 ± 10 
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Figure 2. Population diagram assuming LTE, where Nj is the column den- 
sity of the molecules with the energy state j (Et) and g: the statistical 
weight for Ej. The slope of these lines gives the kinetic temperature of the 
clouds. The upper line represents the 12 CO cloud and the lower line, 13 CO. 



3.3 Large Velocity Gradient Analysis 

The interferometric images of 12 CO (1-0) jKohno et al.|2003| and 
12 CO (2-1) jHsieh et al.||2008f , show that most of CO emission 
comes from the very centre and the starburst nuclear ring of the 
galaxy. In reality, the connection points of the ring with the dust 
lanes of the primary bar present a clear increase of the emission. In 
the following subsections, we analyse the molecular gas physical 
properties of these regions performing a radiative transfer analysis 
of the single dish data described in section|2] 

To improve the estimation of the parameters presented in 
the previous section, we modelled the ratios of the observed lines 
using a non-Local Thermal Equilibrium (Non-LTE) radiative 
transfer code, Radex ( |van der Tak et al.||2007] >, available at 
http : //www . strw . leidenuniv . nl/~moldata/radex .html. 
Radex assumes the Large Velocity Gradient (LVG) approximation, 
and gives the expected line intensity of a molecular cloud with a 
characteristic kinetic temperature T K , molecular Hydrogen density 
n(H 2 ), 12 CO column density Z( 12 CO), and line width v( 12 CO). 

We studied these properties using a x 1 minimization compar- 
ing the observed line ratios and the ratios generated by Radex. We 
built a three-dimensional matrix made of different x 2 values cre- 
ated by changing the three principal input properties in the radia- 
tive transfer code, and varied the inputs within the ranges: 10-150 
K (40 steps) for the kinetic temperature, 10— 10 6 cm 3 (50 steps) for 
the molecular Hydrogen density, and 10 15 — 10 1S cm 2 (30 steps) for 
the I2 CO column density. We fixed the internal 12 CO line width to 
1 km s _1 and the abundance ratio [ 12 CO]/[' 3 CO] to 40 (c.f., 40 ± 10 



found for the centre of active galaxies; Henkel et al. 1993 ; Mauers- 
berger & Henkel|1993||Henkel et al.|1998[|Bayet et al.|2004||Israei| 
2009). We further confirm that changing the abundance ratios be- 



After fixing the abundance ratio and the line width, there are 
only three free parameters (namely temperature, density and col- 
umn density) in the scenario which depict the four observed ratios 
that we are analysing here. However, there are still several possible 
solutions among which we select the best alternative using the com- 
bination of two criteria: the optical depth should be in agreement 
with previous findings (e.g., < 10 as previously found by Petitpas 
|& Wilson[ ;2003 1 and minimizing the combination of the individ- 
ual;^ 2 for every ratio. The errors were estimated using a bootstrap 
method and by looking at the fluctuation of each parameter around 
the^ 2 minimum (see Fig. [5] where the red curves mark the lcr^ 2 
confidence level contours). 



3.4 The Inner 21": One-component Model 

The circumnuclear 21" diameter hosts the active nucleus as well 
as the inner starburst ring with radius « 10". Our single dish ob- 
servations are affected by beam-dilution effect in the measured line 
strength, since the size of the beam is larger than the individual 
clouds, and the dumpiness of the clouds and the large beam size 
produce an unknown filling factor. However, using line ratios with 
an equal beam size, instead of the full line profile, and assuming 
that the emission of every transition comes from the same clouds, 
these problems are avoided (see also Petitpas & Wilson 2003 , |Israel| 
|2009"1 >. 

Applying the x 1 minimization method explained above, we 
generated the confidence levels of the x 1 around the optimal three 
physical parameters: temperature, density and column density (Fig- 
ure We estimated a kinetic temperature of 33+J!? K, a molec- 
ular Hydrogen density of 10 5 "°-»° cm 3 and a 12 CO column den- 
sity of 10 16 ' 7 -o6 cnT 2 , imposing an upper limit for the optical depth 
of r = 4.9 (in agreement with Petitpas & Wilson 2003 ), using an 
abundance ratio [ 12 CO]/[ 13 CO] of 40 and an internal line width of 
1 km s . Upper limit of the optical depth means the biggest value 
estimated by Radex for all transitions. The temperature changes 
by 5 degrees and the molecular Hydrogen density by less that a 
factor 2, when chang ing the [ 12 CO]/[ 13 CO] abundance between 30 
and 50. These estimations are more robust, fully consistent with the 
limits presented by Petitpas & Wilson (2003 1, and present an aver- 
age density larger than 10 4 cnT 3 as expected due from observation 
of HCN(l-O) in the nucleus by jKohno et al.|2003| >. We estimated 
the errors in the magnitudes looking at lcr confidence level and an 
upper limit of the optical depth r < 8. 

Moreover, we obtained a I2 CO column density of 4x 10 16 cm 2 
for a line width of 1 km s 
of 30 km s" 1 



'. If we assume a typical line width 
iKohno et al.|l999|> and a [H 2 ]/[ 12 CO] concentration 



ratio of 10 4 I Wilson et al.|2 009), we obtain a H 2 column density 
of 1.2 x 10 22 cm 4 , which is also in agreement with the estimation 
made using the X C o conversion factor presented in section [3~2j The 
kinetic temperature is in concordance with the value derived using 
the LTE approximation, but its uncertainty of 100 K shows a pos- 
sible contribution of a warmer component. 



tween 30 and 50 introduces only marginal changes in the final re- 
sults. 
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Figure 3. The graphs on the top present the confidence levels of the^ 2 , assuming normally distributed data, around the parameters with minimum^ 2 . The 
minimum value is illustrated by the red star. The red curve marks the posi tion of the lcr^ 2 confidence level. The graphs on the bottom show the upper limit 
optical depth distribution of the lines modelled by Radex, see section [3~3| for details. The red curve marks where the^ 2 confidence level is equal to lcr, and 
the black curve shows the r = 8 contour. The best physical parameters are: T k = 33*\°° K, n(H2)=10 50 -<w> cnT 3 and E(CO)=10 16 ' 7 -<« cnT 2 . 



3.5 Inner 21": Two-component Model 

The central 21"(~1 kpc) of the galaxy hosts a multitude of molec- 
ular clouds with sizes of the order of tens of parsecs and very likely 
exhibiting differing physical conditions. Thus the one-component 
model is not realistic to describe the entire region covered by the 
unresolved central beam of NGC 1097. We take a step further in 
an attempt to describe the typical physical properties of the cold 
interstellar gas using more than one component (cf., Israel 20091. 
We presumed the inner kpc be formed by a cold (30K) and a warm 
(90K) component, and made a similar x 2 minimization study 
explained previously in this section, varying the H 2 densities, 12 CO 
column densities and the fraction between the two components 
T~ = [cold component] /[warm + cold component]. We changed 
the input variables within the ranges: log 10 (n(H2)30K)=[3,6] 
(cnT 3 , 20 steps), log 10 (n(H 2 )90K)=[l-5,4.5] (cnT 3 , 20 steps), 
logio(S( 12 CO))=[15,18] (cm" 2 , 14 steps) and T=[0,1] (8 steps). 
In the selection of the ranges for the densities, we imposed the 
cold component to be denser than the warm component. We fixed 
the abundance ratio [ 12 CO]/[ 13 CO] to 40 and the 12 CO line width 
to 1 km s" 1 , and imposed the same conditions for finding the 
minimum x 2 value as explained previously in this section: a low 
optical depth for the two components and a proper value for the 
X 2 of all individual pairs of line ratios. Our results are shown in 
table 15] with a ratio T = 0.85. It is important to note here that 
the ^ analysis has been carried out using different steps for the 
fitted parameters and different limits for the optical depth and x 1 
of the individual ratios, in order to constrain the confidence of 
the method. The results are correlated within the errors, making 
stronger their confidence. 



Table 3. Kinetic temperature T^, column density n(H2), surface density 
X( 12 CO) and optical depth r from the one- and two-component models. 





1 component 


2 components 

cold warm 


Tk (K) 
log 10 (n(H 2 )) (cm- 3 ) 
log 10 (S( 12 CO)) (cm- 2 ) 
t( 12 CO) 


33 + too 

•"-15 
5.0±>° 

< 4.9 


30 90 

16 7+0.3 17 l ± 0.7 
0.3 2.0 

< 6.6 < 7.3 



4 EMISSION LINES AT DIFFERENT VELOCITY 
REGIONS 

4.1 Dissecting the Inner 21" 

In this section, we use an innovative approach where we "de- 
confuse" the single-dish data to be able to study "resolved" emis- 
sion lines for which no interferometric data have been yet obtained. 
We used the resolved interferometric intensity and kinematic maps 
from |Kohno et al.| {2003} to select five different regions along 
the starburst ring for a deeper analysis of their physical proper- 
ties. Following the combination of two independent criteria, we se- 
lected different knots along the starburst ring to analyse the state 
of the molecular gas. The regions were selected to have a 12 CO(l- 
0) S /N 150 and span a range of velocity greater than one fre- 
quency/velocity channel. These choices have been made to assure 
that the region can be separated from other velocity channels in the 
unresolved single dish data, and have a S IN high enough to return 
reliable measurements for the different transitions . Furthermore, 
all regions have mean line of sight velocities such that they can 
be separated from each other in the velocity channels of the single 
dish data. Two regions are located on the intersection of the bar 
dust lanes with the ring (regions A and B), other two regions at 90 
degrees distance (in the upper and lower side of the ring, regions C 
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Figure 4. The five regions along the inner starburst ring of NGC 1097 with the two bar-ring connection points marked as A and B. All regions are over-plotted 
on the unde rlying 12 CO(1-0) intensity map (Left) . All regions were selected taking into account their distinct velocities and the S/N criteria as described in 
section [4~T| and the full velocity window which contributes to the signal in corresponding channel is marked in the 12 CO(1-0) velocity field (Right). In the 
velocity field, each group of coloured squares indicate pixels with velocity inside the velocity range of each region. 



0.2 



,2 C0(1 -o) 
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Figure 5. The emission lines analysed in this paper (new APEX data, |Kohno| 
|et al.|2003] and |Petitpas~& Wilson 2003 1, showing the contribution to the 
lines of the different regions. 



and D) and a fifth region E in another maximum CO(1-0) emission 
between A and C and close to a connection point (see Fig.[4]l. 

To avoid unnecessary complications regarding the exact shape 
of each star-forming knot, we opted to select a rectangular shape 
along an outer radius of 10.5" and an inner radius between 7-8", 
given their full area fits inside our beam. A narrow velocity range 
is essential for us since we aim to examine the contribution of this 
selected region to the global line, and it is imperative that the dif- 
ferent regions have different velocities to distinguish them in the 
single dish lines. 

Once, the regions were selected and their velocity ranges were 
established, we studied their contribution to the global unresolved 
single dish emission line, see Fig. [5] We took the channels which 
corresponded to the velocity range of each region and generated the 



line ratios using only these channels. We estimated the ratio's errors 
by adding the standard deviation of the channels to the errors listed 
in section[2] With this, we have introduced the error due to a possi- 
ble bad selection of the velocity intervals. Thus, for each region we 
obtained four different line ratios, see table|4] Looking at the ratios, 
it is possible to infer that the physical properties should be simi- 
lar for every region, since the ratios are somewhat similar within 
the errors, but for the case of l2 CO (3-2)/ 13 CO (3-2). To ensure the 
physical properties along the ring, we analysed the ratios for each 
region, following the procedure described in section [33] We real- 
ized the x 1 minimization study described above for every region 
(see section |3~3} , imposing again that the optical depth should not 
exceed 5. The results are shown in the table [4] 

Moreover, we note that the results presented in table|4]are not 
uniquely representative of the physical conditions of the regions 
from the starforming ring only. To estimate the contribution from 
other parts which exhibit velocities which are in the range selected 
for each of the knots along the ring, we calculate the ratio of the 
12 CO(1-0) intensities from the full region inside the velocity range, 
and that from the ring. We find that the contribution of the signal 
from the ring is at best 70% of the pixels inside the chosen velocity 
range (though could be as low as 35%). This fact confirm that single 
dish data cannot be used to rigorously derive the physical proper- 
ties of the ring only, but instead that of regions inside iso-velocity 
contours. Nevertheless, one can assume that the region inside the 
ring (i.e., galactocentric radius < 8.5") is much more homogeneous 
than then ring itself, and the fact that we have derived comparable 
physical properties are indicative of a rather homogeneous ring in 
NGC 1097. 



4.2 The Star Forming Knots Along the Ring 

In star forming rings, stars are expected to form from the accreted 
material because of the shocks and compressions during the ac- 
cretion process (e.g., Elmegreen 1994), and these processes could 
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Table 4. Properties of the different regions inside the central 21" region of NGC 1097. 



A B C D E 



S/N 


151 


152 


151 


152 


150 


Av (km s~') 


41 


46 


57 


19 


26 


12 CO(1-0) 
12 CO(2-l) 


0.87±0.19 


0.83±0.17 


1.3±0.7 


0.77±0.16 


1.0+0.2 


12 CO(3-2) 
12 CO(2-l) 


0.86±0.19 


1.0±0.2 


1.6+0.9 


0.61 ±0.14 


1.0±0.2 


12 CO(2-l) 
13 CO(2-l) 


8.7±1.8 


7.5±1.7 


4±2 


10+1.9 


7±4 


12 CO(3-2) 
13 C0(3-2) 


21±7 


37+18 


15±6 


8.5+1.8 


20±4 


T*(K) 


17±'° 


21+ 10 


124 


21 + 10 
z '-to 


16±j» 


logi (n(H 2 )) (cm- 3 ) 


5 0± 10 


5.4±06 


6.0±1;0 


5.6±°« 


5.0 ± J;° 


log 10 (2:( I2 CO)) (cm- 2 ) 


16.4±0;3 


16.6±JJ;i 


16.3 ± °j 


16.6±°> 


16.4±° 3 



arise randomly by gravitational instabilities without any specific or- 
der along the ring. In the case of NGC 1097, where the ring is pos- 
tulated to be at the location of the inner Lindblad resonance, such 
a scenario would naturally imply that the bulk of the compressive 
forces act on the regions where the ring meets the bar dust lanes, 
and thus the bar may well trigger the star formation in the ring (cf. 
Knapen et al. 2006 Fathi et al. 2008 1 at the connection points and 
the star formation would thereafter successively migrate along the 
ring following the rotation period of a 18 Myr (cf., Bok er et "aT] 



2008). The 12 CO (1-0) and 12 CO (2-1) emission map of the ring 
l |Kohno et al.|2003| | Hsieh et al.|2008| >, present clear peaks of in- 
tensity in the regions around the dust lanes connection points with 
the ring (A, B and E in Fig. ffl, and hence, may indicate stronger 
star formation on those spots. This could imply an increase in the 
molecular gas density with respect the rest of the ring. 

Our improved LVG analysis using five different regions at dif- 
ferent velocities, does not show significant differences between the 
molecular gas densities, temperatures or column densities, see ta- 
ble|4] It is possible that the measured similarity between the derived 
molecular gas densities, could be due to a bad estimation of the con- 
tribution of the ring to the observed line. However, the region inside 
the ring (radius < 8.5") shows relatively small 12 CO(1-0) intensity 
variation as well as small variations in Ha intensity as shown in 
high resolution Hubble Space Telescope images presented in Fathi 
|et al. | p006 1. It is thus likely that the star forming gas inside the ring 
is rather homogeneously distributed, hence, its contribution to our 
derived physical properties along the different regions described in 
section |4~T| is similar. Consequently, our derived values hint that 
the conditions of the molecular gas at the different regions along 
the ring are comparable. 

These results are consistent by the recent far infrared and sub- 
mm studies, from the Herschel Space Observatory, of the ring in 
NGC 1097 which showed no azimuthal age gradient (Sa ndstrom| 
|et al.|2010) , accompanied by an increase of the ionized gas density 
in over the full ring associated with star formation activity (Beirao 
|et al.|20i0) . Moreover, we compare the regions selected here with 
the resolved Very Large Array radio continuum maps (Bec k et aT] 
|2005| their Figure 14) and find that all star forming knots exhibit 
radio emission, indicative of ongoing star formation at comparable 
intensities. 

The strong 12 CO (1-0) and 12 CO (2-1) emission observed 



with the interferometric data JKohno et al.| [2003 and |Hsieh| 
|et al.| [2008^ on the connecting points may be due to partic- 
ularly strong shocks, as predicted by the peak of the ratio 
[OI]63pm/[CII]158yum ( |Beirao et aT1|2010^ , and not necessarily 
due to triggered strong star formation. 



5 MASS INFLOW RATE ONTO THE SMBH 

The fuelling of the active nucleus of NGC 1097 has been subject 
to much debate recently. Modelling the spectral energy distribu- 
tion, [Nemmen^enir] ((2006 ) estimated a inflow mass rate of m = 
6.4 x 10 3 M E[M , with the Eddington accretion rate of M Edd = 2.7 



M yr 1 , for a 1 .2x 10 s M Q black hole. Recently, [van de Ven & Fathi 



( |20 10| > used two-dimensional spectra in the optical to kinematically 
derive the inflow rate from ionized gas. These authors measured 
the inflow velocity from a kinematic analysis and the electron den- 
sity from [Sn] observations, consistent for the most inner part with 
Nemmen et al. ( 2006 I. Davies et al. ( 2009 I estimated an upper limit 
1.2 M yr -1 . These studies are based on simplifications of the in- 
flow mechanism along the nuclear spiral arms and also on indirect 
measurements of the inflowing material by assuming that the hot 
ionized gas traces the underlying cold molecular gas. 

We estimated the mass inflow rate using our direct measure- 
ment of the molecular gas density, p gas , obtained for the central kpc 
of NGC 1097 (see section [3~4} . We combined it with the inflow ve- 
locity estimation from | van de Ven &"F athi (2010), and their mass 
inflow rate analysis. 

, h 1 

M = m v inflow Ap nR- - — , (1) 
R Am 

where m is the number of arms, v m ^ „ the inflow velocity, Ap gas 
the over-density on the arms, R the galactocentric radius and h the 
disc scale height. We assumed a nuclear spiral arms over-density 
of the 10%, i.e., Ap gas = 0.1p gas (e.g. |Englmaier & Shlosman|2000| 
|Maciej ewski 2004), and for comparison also measured the mass 
inflow rate using the electron over-density measured by |van de Ven| 
|& Fathi| ( |2010} . We further compare the results estimating the scale 
height h taking into account the thickness of the disc, and used the 
criterion of marginal stability geff = 1 (Ro meo|1994| ), where 

geff ~ 2d + 20; t= x^L i Q=-t^- ( 2 ) 
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Figure 6. Mass inflow rate of the inner 400 parsec radius of NGC 1097. 
Solid line mass inflow rate estimated by |van de Ven & FalE|(2l)To). Green 
line m ass inflow rate estimated by |NemmeiTet*^l . (20061. Nemm en et al.| 
{2006} assumed a mass inflow rate of M = M (r/r tr ) p , where p = 0.8, 
r lr = 225 and M D = 6.4 X 10~ 3 MEdd- Dotted line our mass inflow rate 
estimation from the total molecular gas and using the measured of the elec- 
tron over-density of the arms by | van de Ven & F athi 1 2010 ». Dashed line 
our mass inflow rate estimation from the total molecular gas and using an 
over-density of the arms of the 10%. Red dashed line our mass inflow rate 
estimation from the total molecular gas and using an over-density of the 
arms of the 10% and using a scale height from Romeo ( 19941. Red dotted 
line our mass inflow rate estimation from the total molecular gas and using 
the measured of the electron over-density of the arms by jvan de Ven & Fathi 
(20i0) and using a scale height from Romeo ( 1994}. The mass inflow rate 
estimated using the molecular gas density and the electron over-density for 
a thin disc at around 70 pc is 0. 10 M Q /yr. 



k is the epicyclic frequency, c s the sound speed (we use c s = 
10 km s~') and E gas « Pg^/h the gas column density. We did 
not find a significant difference between this assumption and the 
marginally stable disc, and we attribute this to the fact that we 
are measuring cold molecular gas which is thinner than its ionized 
counterpart, i.e., the cold gas component has h/R < 1. We obtained 
several measurements for the mass inflow rate onto the SMBH in 
NGC 1097 as shown in Figure [6] 

Our measurement using the thick disc description is 0.11 
M /yr at 1"= 70 pc distance from the SMBH, which is the ra- 
dius to which the line of slight velocities were reliably measured 
by |van de Ven & Fathi (2010). This value is still in agreement with 
the accretion estimation for the transition between Seyfert 1 and 
Low-ionization nuclear emission-line region galaxies jHo]|2005|, 
however, puts NGC 1097 on the more efficient end of the distri- 
bution of the transition objects. Our inflow rate at 70 pc is within 
the upper limit derived by Davies et al. ( 20 09) and is one order of 
magnitude greater than the value derived by |van de Ven & Fathi| 
l |2010} .The difference with the latter work is driven by the fact that 
while our cold disc remains thin (h/R < 1), the ionized disc of |van| 
|de Ven &Fa thi (2010) thickens and to avoid, this these authors ap- 
plied an upper limit for the h/R, leading to a lowering of the mass 
inflow rate (c.f. equation[TJ. 



6 X co CONVERSION FACTOR 

With the analysis of marginal stability of the disc presented in sec- 
tion]^] we have been able to estimate from the kinematics analysed 
by |van de Ven & Fathi (2010), the column density of the gas along 



the nucleus using |Toomre| ( [T964) (cf. |van de V en & Fat hi|2010| ) and 
|Romeo| ( [l994[ > criteria. This method for estimating the surface den- 
sity is completely independent of the CO emission in the galaxy. 
Assuming Q=l (or Q e // = 1) and given cv,=10 km s~ ] as well as k 
(from the "kinematic model" of van de Ven & Fathi (2010)), yields 
a prediction for E gfl! ; then assuming S S(1I = Z// 2 , the ratio with the 
CO intensity yields the X C o factor. Hence, using this calculated 
column density for different radii, we could in principle study the 
Xco conversion factor for a given I2 CO (1-0) intensity measured 
with the corresponding beam size. We translate this surface col- 
umn density to molecular Hydrogen column density, dividing by a 
factor 1.36 corresponding to the Helium contribution to the gas. 

We have the l2 CO (1-0) integrated velocity intensity for a 
beam size of 15" I Ko hno et al.|2003l and 21" (see section^, which 
correspond to radius of 7.5"and 10.5" , respectively. Hence, we 
calculated the average surface density determined by the marginal 
stability criteria, over an area similar to the beam size. Combin- 
ing the surface densities and 12 CO (1-0) intensities for equal ar- 
eas, we estimated the X co conversion factor. The results were 
X co = (2.8 ± 0.5) x 10 20 cm" 2 (Kkms" 1 )- 1 at radius of 10.5" and 
Xco = (5.0 ± 0.5) x 10 20 cnr 2 (Kkms" 1 )- 1 at radius of 7.5". For 
the estimation of the X C o error, we took into account the contribu- 
tion to the uncertainty of the epicyclic frequency k (uncertainty of 
7%), and of the 12 CO (1-0) integrated velocity intensity for each 
beam (3% for a beam size of 15" and 10% for 21"). Moreover, for 
these calculations we assumed a sound speed equal to 10 km s -1 , 
which is a lower limit for the sound speed consider by |Davies et al.| 
( 2009 1, 20-40 km s -1 , what makes our X C o > also a lower limit. 

We obtained a X C o conversion factor comparable to that pre- 



sented in section [3T2] Moreover, the surface density estimated by 
the marginal stability analysis is (2.4±0.2)x 10 22 cirr 2 , which again 



is in agreement with our estimations in section |32| and section |3T4| 
While our derived Xco conversion factor compares well with the 
typical values assumed in the literature, it is in disagreement with 
other estimations of this factor in centre of galaxies (e.g., [Strong] 
|et al.|2 004 Israel 2009 and references therein), which are one or- 
der of magnitude lower. This disagreement could be due to the fact 
that we have not accounted for the ionized and atomic gas column 
density in the total gas column density estimation. 



7 CONCLUSIONS 

We have analysed single dish observations of five different tran- 
sitions of the CO molecule in the central 21" (kpc) region of 
NGC 1097. We have used X co conversion factor, LTE and LVG 
approximation to parametrize the cold ( « 30 K) phase of the in- 
terstellar gas from which we have identified and characterized the 
typical physical properties inside this region. All methods have re- 
sulted in comparable values, and in agreement with other estima- 
tions made by |Petitpas & Wilson| ( [2003] ) and |Hsieh efaT] ( [2008) , 
who used a different set of line ratios and a different LVG approx- 
imation code. We obtained a typical kinetic temperature of about 
33 K, a molecular Hydrogen density of 4.9x 10 3 M G pc~ 3 and a CO 
column density of 1.2 x 10~ 2 M Q pc~ 2 . 

We have further applied a scheme which involves the assump- 
tion that circumnuclear kpc region of NGC 1097 consists of two 
dominant components, a cold and a warm («; 90 K) component (cf. 
Israel 2009 ). The two component analysis yields a 85% cold/total 
ratio, where the cold gas is denser than the warm component. This 
ratio and the physical properties of the cold and warm gas are sim- 
ilar to the values found in NGC 1068 ( |Israel|2009> . 
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To gauge the robustness of these results, we dissect the spa- 
tially unresolved observations over the central kpc of NGC 1097 
in velocity space. We apply the LVG analysis on line ratios de- 
rived from five carefully selected velocity channels and find that the 
physical properties remain unchanged to within the errors. These 
five regions were selected based on the interferometric 12 CO(1-0) 
intensity and velocity map, and all contain relatively strong star 
forming knots of the star burst ring. The contributing signal from 
the star forming knots is at least of S/N = 150 and the contribu- 
tion of each knot to the full 12 CO(1-0) within the corresponding 
velocity bin varies between 35 and 70%. Based on the interfero- 
metric data and multi-wavelength images in the literature we can 
assume that the cold interstellar gas in the region interior to the ring 
is fairly homogeneous, thus find it likely that the physical parame- 
ters derived from these regions relect the inherent physical proper- 
ties along the starburst ring of NGC 1097. Using the Schmidt law, 
the comparable gas physical properties imply similar star formation 
rates along the ring. Although we do not make any strong claim 
about this fact, it is encouraging that they are in agreement with 
the recent data from the Herschel Space Observatory, supporting 
an star formation scenario described by Elmegreen ( 19941. 

The rotation curve presented in |van de Ven & Fathi] j2010) 
which reaches approximately 400 km s~' at 8" radius, yields the 
orbital time of 18 Myr at this radius. In the simplified case, where 
the ring is fed by clumps of material from the bar dust lanes, fol- 
lowed by episodic star formation in the ring connection points with 
the bar dust lanes, the star forming complexes could migrate along 
the ring in a very short time. In such a case the star forming knots 
along the ring are expected to exhibit similar properties, since they 
"mix" very efficiently. This is further supported by the fact that the 
knots coincide with knots of enhanced radio continuum emission 
l |Beck et al.|2005) . 

Making use of our new derivations of the cold gas densities, 
we have been able to revise the mass inflow rate onto the SMBH 
in NGC 1097. We derived an inflow rate of 1.1 M G /yr for a thick 
disc, and conclude that accounting for the total interstellar medium 
and applying a careful contribution of the disc thickness and cor- 
responding stability criterion increases the inflow rate by a factor 
10. The critical value 0.01 Mem for the transition between LINER 
and Syl galaxies was found by Ho (2005 1 with a distribution over 
three orders of magnitudes. Our new measurement of the accre- 
tion rate onto the SMBH of NGC 1097 places this galaxy at the 
higher end of this distribution, and combining our derived accre- 
tion rate with the updated Eddington ratio presented in |Eracleous| 
[etlTlpOTol l, NGC 1097 is thus placed in the regime between the 
advcction-dominatcd mass accretion and accretion via a thin disc 
structure (e.g., Narayan et al.|1998| >. 

We have moreover calculated the Xco conversion factor for 
the centre of NGC 1097 using an independent estimation of the 
surface density to the CO emission. We used the same analysis 
described in |van de Ven & Fath i (20101 from the kinematics of 
NGC 1097. We obtained X co = (2.8±0.5)x 10 20 cm~ 2 (Kkms-'r 1 
at radius 10.5" andX co = (5.0±0.5)x 10 20 cnT 2 (Kkms-'r 1 at ra- 
dius 7.5". The surface density given by this "kinematical method" 
and the X co conversion factor are in agreement with our estima- 
tions, see section [3~2] and section [3T4| However, these values found 
at the centre of NGC 1097 are in disagreement with values found 
for centres of other active galaxies which are one order of magni- 
tude lower (e.g. |Strong et al.|2004||Israel|2009} . 

With the approach and analysis described in this paper we 
have demonstrated that important physical properties can be de- 
rived to a resolution beyond the single dish resolution element. 



To obtain a more accurate description of the molecular gas phys- 
ical properties of the nuclear region of NGC 1097, observations of 
higher transitions, e.g., CO (4-3), with the Atacama Large Millime- 
tre Array (ALMA) will help overcome most of the limitations that 
we have faced in the current work. The ALMA data will (i) give 
us the complete description of the properties of star forming knots 
in the ring, (ii) will allow us to calculate the exact mass inflow rate 
for different radii by resolving the density for different radii, and 
(Hi) will resolve the unknown nature of the Xco conversion factor 
in the circumnuclear kpc of this prototype LINER/Syl Galaxy. 
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